Evaluation of the MODIS retrievals of dust aerosol over the ocean during PRIDE 
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Abstract: 


The Puerto Rico Dust Experiment (PRIDE) took place in Roosevelt Roads, Puerto Rico 
from June 26 to July 24, 2000 to study the radiative and physical properties of African 
dust aerosol transported into the region. PRIDE had the unique distinction of being the 
first major field experiment to allow direct comparison of aerosol retrievals from the 
MODerate Imaging Spectro-radiometer (MODIS) with sunphotometer and in-situ aerosol 
measurements. Over the ocean, the MODIS algorithm retrieves aerosol optical depth 
(AOD) as well as information about the aerosols’ size distribution. During PRIDE, 
MODIS derived AODs in the red wavelengths (0.66 pm) compare closely with AODs 
measured from sunphotometers, but, are too large at blue and green wavelengths (0.47 
and 0.55 pm) and too small in the infrared (0.87 pm). This discrepancy of spectral slope 
results in particle size distributions retrieved by MODIS that are small compared to in- 
situ measurements, and smaller still when compared to sunphotometer sky radiance 
inversions. The differences in size distributions are, at least in part, associated with 
MODIS’ simplification of dust as spherical particles. Analysis of this PRIDE data set is a 
first step towards derivation of realistic non-spherical models for future MODIS 
retrievals. 


2 


1. Introduction 


Mineral dust aerosols are produced mainly by wind erosion of desert soils, and are a 
significant component of tropospheric aerosols [Prospero, 1996; Chiapello, et al. 1999]. 
These aerosols are lifted by the wind, raised to high altitudes by convection, and may be 
transported over long distances from their sources [Ginoux et al., 2001 , Li- Jones and 
Prospero, 1998; Formenti et al., 2001; Smirnov et al., 2000b]. They influence the optical 
properties of the Earth’s atmosphere and climate through the scattering and absorption of 
sunlight [Tanre et al., 2001], thereby in turn influencing local and global atmospheric 
dynamics [Alpert et al., 1998; Miller and Tegen, 1998]. As for indirect climate forcings, 
dust aerosol may influence photochemical processes [Dickerson et al., 1997], contribute 
to cloud condensation nuclei [Levin and Ganor, 1996], or act to suppress precipitation 
[Rosenfeld et al., 2001]. Dust is deposited into the ocean [Gao et al., 2001], and is related 
to the biological productivity of a basin. Mineral dust in large quantities affects visibility, 
and may have adverse influences on human health [Prospero et al., 1999], and upon 
animal populations [Stallard, 2001]. Recent studies, mentioned by Sokolik et al. [2001], 
suggest that some change in dust production may be caused by anthropogenic activities 
[Tegen et al., 1996; Tegen and Fung, 1995], though satellite data analysis indicates 
dominance of sources in scarcely populated regions [Prospero et al., 2002]. 

A major source region for mineral dust is the combined Sahara and Sahel areas of North 
Africa [ Prospero , 1996; Moulin et al., 1997; Johansen et al., 2000]. This source is active 
nearly all year, with plumes flowing across the Atlantic towards the Caribbean and the 
Americas, especially during the summer months [Prospero, 1996; Higurashi et al., 

2000]. These plumes have been well documented by satellite sensors such as the Total 
Ozone Mapping Spectrometer (TOMS) [Herman et al., 1997; Chiapello et al., 1999] and 
the Advanced Very High Resolution Radiometer (AVHRR) [Husar et al., 1997; 
Higurashi et al., 2000]. Ground based in-situ instruments and sunphotometers have also 
observed dust in the Caribbean, at sites such as Barbados [Prospero, 1996; Smirnov et al., 
2000b]. 
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For some dusty regions of the world, satellite observations have been linked to 
sunphotometer data in the works of Higurashi et al. [2000] and Livinsgston et al. [2000] 
(using AVHRR data), Tanre et al. [2001] (using Thematic Mapper data.), and Moulin et 
al. [1997] (using Meteosat data). Specifically, in the Caribbean, satellite retrievals have 
been compared with in-situ measurements by works such as Husar et al. [1997] and 
Chiapello et al. [1999]. These studies have shown that satellite retrievals are a promising 
method for identifying dust and retrieving its properties. 

As part of the NASA’s Earth Observing System (EOS), the Terra satellite was launched 
in December 1999. Its flagship instrument, the MODerate resolution Imaging 
Spectrometer (MODIS) [ Salmonson et al., 1989; King et al., 1992] measures the ambient 
aerosol column over oceans at seven wavelengths in the solar spectrum, with high 
accuracy and high resolution over a variety of time scales. Aerosol algorithms for 
MODIS over ocean [Tanre et al., 1997] and over land [Kaufman et al., 1997] have been 
validated under a variety of conditions [Remer et al., 2002; Chu et al., 2002], but no 
evaluation has been made specifically for dust over the ocean. 

Retrieving dust aerosol properties from satellite may prove to be particularly challenging, 
mainly due to the non-spherical shape of dust particles. Microscopic analyses of dust 
particles show that they are irregular in shape rather than spherical [Koren et al., 2001]. 
However, the dust optical properties for MODIS, similarly to previous satellite missions, 
are modeled by assuming that the particles are spherical [Tanre et al., 1997]. 

Investigators such as Mischenko et al. [1997] showed that non-sphericity may have large 
effects on the scattering optical properties of the aerosol, especially at large scattering 
angles (greater than 120°) that would be seen by MODIS. 

An opportunity for dust validation arose with the Puerto Rico Dust Experiment (PRIDE) 
[Reid, 2000; Reid et al, 2002]. Held from June 26 to July 24, 2000, PRIDE was operated 
from Roosevelt Roads, Puerto Rico, and was designed to study the African dust aerosol 
transported into Puerto Rico. For MODIS, PRIDE had the unique distinction of being the 
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first major field experiment to allow direct comparison of its aerosol retrievals to field 
measurements. 

In this study, we report on the evaluation of the MODIS aerosol retrievals from MODIS 
during PRIDE, using sunphotometer and in-situ observations. In section 2, we outline the 
theoretical and operational use of the MODIS aerosol over ocean algorithm. Section 3 
describes each instrument used to validate the MODIS retrievals, and the data taken 
during PRIDE. In section 4, we show the comparisons of MODIS and validation data, 
and in section 5, we discuss how these comparisons should be used as a basis for new 
science on dust aerosol. 


2. MODIS Aerosol Retrieval 

The MODerate Imaging Spectro-radiometer (MODIS) is the flagship instniment aboard 
the Earth Observing System (EOS) satellites [Salmonson etal., 1989; King etal., 1992]. 
MODIS performs measurements at 36 channels in the solar and infrared regions (0.415 to 
14.235 pm), with resolutions of 250 m, 500 m, or 1 km, depending on the wavelength. 

At a nominal altitude of about 700 km, MODIS observes a swath about 2300 km wide. 
The first MODIS instrument was launched with the Terra satellite in December, 1999, 
which has a sun-synchronous orbit that passes southward over the equator at 10:30 AM 
local (solar) time. Aboard Terra, MODIS provides nearly global coverage each day. 
MODIS ’s wide spectral range and fine spatial resolution, coupled with its broad swath 
over the Earth’s surface, make it suitable for monitoring events on short term local time 
and spatial scales, as well as for global and long term scales. Previous papers by 
Kaufman et al. [1997] and Tanr4 et al. [1997] discuss using MODIS to retrieve aerosol 
properties over the land and ocean, respectively. 

2.1. Theoretical Description of the Retrieval Algorithm 
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Tanre et al. [1997] details the strategy for using MODIS to retrieve aerosol properties 
over the ocean. Observed top of the atmosphere (TOA) reflectances at six wavelengths 
(0.55, 0.66, 0.87, 1.24, 1.64 and 2.13 pm) are compared with a lookup-table of pre- 
computed reflectance for an array of angles, size distributions and optical depths. The 
modeled reflectance with the smallest difference from the observed reflectance is 
retrieved from the look up table. This best fit reflectance is associated with a 
corresponding set of aerosol properties, which are considered to be the retrieved products. 

The reflectances in the lookup table are computed from aerosol models that represent the 
aerosol properties of a vertical column. Currently, nine tropospheric aerosol models are 
used, including four “fine” mode models (accumulation mode: dominated by chemical 
and combustion processes) and five “coarse” mode models (dominated by maritime 
particles and dust). The optical models are described in Tables la and lb. The current 
lookup tables were updated from those described by Tanre et al. [1997], the main 
difference being the use of new “dust-like” particle models (large modes 4 and 5). Their 
scattering and absorption properties were derived from a combination of AERONET data 
and LANDS AT Thematic Mapper (TM) image analysis [Tanre et al., 2001; Kaufman et 
al, 2001], 

For each model, the modeled satellite signal is assumed to be a combination of radiation 
from the atmosphere and reflection from the surface. The atmospheric calculation 
accounts for multiple scattering by molecules and the aerosol, as well as reflection of the 
atmosphere by the sea surface. The ocean surface calculation includes three 
contributions: the Fresnel (“sun glitter”) reflection off the surface waves, reflection by 
whitecaps and foam and Lambertian reflectance coming from underwater scattering 
(sediments, chlorophyll, etc). The surface wind speed (for sunglitter and foam 
calculations) is assumed fixed at 6.0 m/s. Zero water leaving radiance is assumed at all 
wavelengths, except for at 0.55 pm, where a reflectance of 0.005 is used. 

Using the radiative transfer code developed by Ahmad and Fraser [1981], spectral 
reflectances were computed for each of the nine aerosol models. Five values of aerosol 
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columnar optical depth, x a , (total aerosol loading) at 0.55 pm are considered for each 
mode, ranging from a pure molecular (Rayleigh) atmosphere (x a = 0.0) to a highly turbid 
atmosphere (x a = 2.0), with intermediate values of 0.2, 0.5 and 1.0. For each model and 
aerosol optical depth at 0.55 pm, the associated optical depths are stored for the other five 
wavelengths, plus an additional wavelength in the blue (at 0.47 pm). Computations are 
performed for combinations of 9 solar zenith angles, 16 satellite zenith angles and 16 
relative sun/satellite azimuth angles (2304 total combinations). 

To perform the aerosol retrieval, we use the method discussed by Tanre et al. [1997]. 

The multiple scattering radiance from two lognormal modes (we assume one small, one 
large) can be approximated by the weighted average of the two modes, calculated for the 
same optical thickness [ Gordon et al., 1997]. Let us assume that the total reflectance 
measured at a wavelength channel X, p x m , (superscript “m” denotes “measured”) at the 
satellite level is: 

p:=r)p[+( '-r])p[, ( 1 ) 

where p x s and p>' are the reflectances of the small (s) and large (1) modes, respectively, 
and r) is the ratio of the reflectance contribution from the small mode to the combined 
modes. These reflectances include the contributions from the surface and rayleigh 
scattering. 

Given one of the 20 combinations of small and large modes, we compute the expected 
reflectance by interpolating to the exact sun/satellite geometry, for each of the aerosol 
optical depths in the table. Starting at 0.87 pm, we derive the exact AOD for the 
combination by comparing to the observed reflectance, and then using this to derive the 
optical thickness at the other wavelengths. From the two modes, the small mode ratio, r|, 
is computed, which is the ratio of the optical depth contributed by the small mode to the 
total optical depth. 

T| = T s / X (2) 
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If we denote the calculated reflectance at channel X by p x csl , (where si represents the 
combination of small and large modes) and the measured reflectance by p x m , then the 
relative error between the two, e x sl , is given by: 


e , = 


p: -p? 

pr'+o.oi 


(3) 


where the constant 0.01 is a small residual used to prevent division by zero. The total 
relative error for all wavelengths, e sl , is given by: 


si 


e 



(4) 


where the relative error from the blue (0.47 pm) is not included in the calculation at this 
time, due to potentially large uncertainties in oceanic pigment concentrations . Whichever 
ratio T), of the combination of modes, si, gives the smallest total relative error is 
considered the “best” solution. For operational purposes, we often use the average for the 
three combinations with the smallest error. Going back to the lookup tables, we retrieve 
the number size distribution for each mode, and weight them to compute the effective 
radius R^of the particle population, defined as: 

tor 

(5) 


R „- 


r 
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r 
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where r is the radius. 


2.2. Retrieval Algorithm in Operation 

Proper preparation of the MODIS input is necessary for valid aerosol retrievals over 
ocean. MODIS measures its 36 channels at three different resolutions. Using a geo- 
location file, MODIS data are separated into land and ocean. For the ocean aerosol 
retrieval, reflectance inputs include: two wavelengths (0.66 and 0.87 pm) at 250 meter 
resolution, and five wavelengths (0.47, 0.55, 1.24, 1.64 and 2.13 pm) at 500 meters. 
These reflectances are separated into 10 km x 10 km boxes (i.e. 20 x 20 pixel boxes for 
the 500 m channels, and 40 x 40 pixel boxes for the 250 m channels). Using a 
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combination of reflectance thresholds and variability tests [Martins et al., 2002], data are 
quality screened and a cloud mask is produced. 

Using the reflectance at 0.87 pm, we remove the brightest 25% and darkest 25% of the 
remaining (non-screened) pixels in each 10 km x 10 km box, in order to minimize 
contamination arising from inhomogeneous surfaces and subpixel cloud and cloud 
shadow features. Averages and standard deviations of the remaining pixels are stored. 

Before doing the actual aerosol retrieval, we perform a number of tests on the reflectance 
data; for example, ensuring that all reflectance values are at least the minimum rayleigh 
value. Additionally, we check the glint angle, defined as: 

0gii m = cos _1 ((cos0 s cos0 v )-(sin0 s sin0 v cos(J))), (6) 

where 0 S , 0 V , and 4> are the solar zenith, the satellite zenith and the relative azimuth angles 
(between the sun and satellite), respectively. To avoid glint contamination, we do not 
retrieve aerosol properties on boxes within 40 degrees of the specular reflection angle. 

Remer et al. [2002] has shown that, as compared with oceanic AERONET sunphotometer 
data, the ocean algorithm gives optical depth accuracy of Ax - ± 0.03x ± 0.05. For the 
effective radius, the accuracy was calculated to be ±0.10 ftm. However, all of their.' 
validation points were of non-dust aerosol. 

2.3. PRIDE Data 


Table 2 lists the date and time of each MODIS overpass over Roosevelt Roads, Puerto 
Rico during PRIDE. Included are data for the solar and satellite zenith angles and the 
relative azimuth angle between the two (all measured from the surface). The next 
columns define the path scattering angle, and the computed glint angle (from equation 6). 
The last column denotes whether Roosevelt Roads is within the defined glint masked 
region (40 degrees from the specular angle). Note that MODIS does not operationally 
retrieve aerosol properties over the ocean for “Glint” overpasses. 
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Some of the higher optical depth conditions during PRIDE occurred on July 21, 2000 
(Julian day 203), when the optical depth at 0.55//m was about 0.5. Figure 1 shows the 
MODIS “true-color” image, produced by combining reflectance data from the MODIS 
red, green and blue channels (0.66, 0.55 and 0.47 pm). In the image, Puerto Rico (at 
.1 8°N, 66°W) is located about one-third of the way between the left and right sides, . 
approximately centered vertically. Notice the visible glint centered near the top of the 
image. 

Figure 2 is the corresponding MODIS aerosol retrieval at 0.55 pm for the image in Figure 
1 , combined for land and ocean. The potential “glint” region (40 degrees) is masked out 
in gray, along with areas removed by the cloud mask. This 40° glint mask is purposely 
conservative, and encompasses a much larger area than the area anticipated by the visible 
glint in Fig 1 . Also, some of the glint may not be visible in Fig l because it is obscured 
by heavy dust. At the time this data was acquired, the main pulse of the high optical 
depth was near or slightly west of Puerto Rico. To the east of the “glint” mask, we see 
evidence of the moderate dust plume that hit Puerto Rico a few days later. . . 


3. Validation Data Sets 

For validating MODIS data, we employed data from multiple platforms. For direct 
comparison of optical depth, we used three different types of sunphotometers. For 
analysis of MODIS size distribution retrievals, we used both AERONET sunphotometer 
almucantars and ground-based in-situ retrievals. As a schematic illustration, figure 3 
shows where instruments were deployed in relation to the MODIS satellite track on 4 
July. 

3.1. Sunphotometer data 
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For the validation of MODIS optical depth, we use data collected from three types of 
sunphotometer instruments. Two automatic sunphotometers were provided by the 
AERONET program [Holben et al., 1998] for the duration of PRIDE. Two 
MICROTOPS II handheld sunphotometers [ Morys et al., 2001] were deployed at various 
times and locations, including during ship cruises. Finally, the 6 channel Ames Airborne 
Tracking Sunphotometer (AATS-6) [Matsumoto et al., 1987] was flown aboard the 
Navajo aircraft on numerous flights [Livingston et al., this issue]. On different days, 
measurements from some or all of the sunphotometers coincided with MODIS overpass. 
Spectral ranges are 0.44 - 1.02pm for AERONET, 0.34 - 0.87pm for MICROTOPS and 
0.38 - 1.02pm for the AATS-6. 

3.1.1. AERONET Sun/Sky-photometers 

The AERONET network is a global network of sun/sky autonomous radiometers. 
Description of the instruments and data acquisition is given by Holben et al. [1998] . In 
“sun” mode, the instrument automatically tracks the sun, retrieving optical depth from 
measurements of solar extinction. In “sky” mode, the instrument measures radiance in 
the sun’s almucantar (same zenith angle, varying azimuth), which are later used to 
retrieve aerosol size distribution and other parameters. 

Two such instruments were deployed for the duration of the PRIDE experiment, one 
along the east coast of Puerto Rico at Roosevelt Roads (Latitude = 18.20 C N, Longitude - 
65.60°W) the other along the south coast in La Paguera (Latitude = 17.97°N, Longitude = 
67.05°W). These instruments (especially Roosevelt Roads) were expected to measure 
dust directly transported over the ocean, with little contamination from other sources. 

During PRIDE, both instruments performed sky and sun measurements in at least four 
spectral bands, 0.44pm, 0.67pm, 0.87pm and 1.02pm, from which the aerosol optical 
depth and the aerosol size distribution were derived. Direct optical depth measurements 
were taken approximately every 15 minutes, while sky radiance data were observed every 
hour. Figure 4 presents a time series of PRIDE optical depth at both Roosevelt Roads 
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and La Paguera. Optical depth calibration coefficients are based on inter-comparison with 
a reference instrument (usually one calibrated at the pristine mountain top of Mauna Loa, 
in Hilo, Hawaii), whereas radiance measurements are calibrated using an integrating 
sphere at NASA-Goddard Space Flight Center. Expected error is less than 0.01 for the 
optical depth, and less than 5% for the sky radiance [Holben et al., 1998]. 

Direct optical depth measurements used in this study, have been quality checked and 
cloud screened using the method of Smirnov et al. [2000a]. Quality controlled radiance 
measurements were turned into retrievals of aerosol size distribution by the method of 
Dubovik and King [2000]. The third column of Table 3 lists almucantars taken within two 
hours of MODIS overpass (later used for size distribution analysis). 

3.1.2. Handheld Sunphotometers 

Two handheld sunphotometers (MICROTOPS II) were deployed for PRIDE. These 
instruments, manufactured by Solar Light Co, Inc (Philadelphia, PA USA), weigh less 
than a kilogram and measure 10 x 20 x 4.3 cm [Morys et al., 2001]. Both of our 
instruments were identical, intended to measure AOD at four channels (0.340, 0.440, 
0.675 and 0.870 pm) and the precipitable water column using 0.936 pm. 

Calibration for these instruments was done from inter-comparison with a reference 
sunphotometer [Ichoku, et al., 2002b], Usually this instrument was a reference 
AERONET instrument located at Goddard Space Flight Center (regularly calibrated by 
performing Langley plot analyses atop the pristine Mauna Loa Observatory in Hilo, 
Hawaii). Calibration error for MICROTOPS optical depth is no more than 0.02 [Ichoku, 
et al., 2002b], 

Tests show that two main sources of error are improper pointing at the sun, and improper 
cloud screening decisions [Ichoku et al., 2002b], In practice, users of these MICROTOPS 
II sunphotometers undergo proper orientation before use, in order to minimize these 
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human errors. Even so, we require three or more scans in quick succession, per 
observation, and retain only those triplets that are sufficiently constant in value. 

These instruments are highly portable and can be deployed at locations where the 
logistics of installing other types of sunphotometers may be impossible. During one three 
day period, intensive observations were taken aboard the University of Puerto Rico’s RV 
Chapman, at various ocean locations east of Puerto Rico. We learned that taking 
measurements aboard a rolling ship causes unique problems for sun-pointing, but we 
believe that taking multiple scans per observation enabled quality control. The fourth 
column in Table 3 lists the dates and location of calibrated MICROTOPS measurements, 
corresponding to MODIS overpass. Observations labeled “Roosy Roads” were taken at 
several locations, but always within a few kilometers of the AERONET sunphotometer at 
Roosevelt Roads. 

3.1.3. Ames Airborne Tracking Sunphotometer 

The Ames Airborne auto-Tracking Sunphotometer automatically tracks the sun and 
measures the transmitted the solar beam, retrieving the overlying columnar optical depth 
[Russell et al., 1993; Matsumoto et al., 1987]. During PRIDE, the 6-channel version was 
mounted in the Navajo aircraft and was set to observe AOD and column water vapor. It 
measured direct solar beam transmission in six spectral channels (380.1, 450.9, 525.7, 
864.5, 941.9, and 1021.3 nm, with filter bandwidths of 5.0-5. 8 nm) [Livingston et al., this 
issue]. Data were digitized and recorded every 3 seconds. 

The AATS-6 was calibrated before and after PRIDE, by taking it to the pristine 
conditions of the Mauna Loa Observatory (in Hilo, HI), and performing Langley plot 
analyses [Russell et al., 1993]. AATS-6 data quality during PRIDE are discussed in 
Livingston et al., [this issue]. 

During PRIDE, the Navajo flight paths were designed in part so that the AATS-6 could 
measure the vertical distribution of AOD. However, most flights included a level track 
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within 100 meters of the ocean surface, between 5 minutes before and 5 minutes after 
MODIS overpass. By assuming that only a small fraction of the total optical depth was 
located below the aircraft, we used the statistics from the AATS-6 during over-pass time 
to compare with MODIS data. The fifth column of Table 3 lists the AATS-6 flights 
during PRIDE, where the latitude and longitude values are the midpoint of the near- 
surface flight track. 

3.2. In-situ Measurements 

During PRIDE, ground based aerosol sampling and in-situ optical measurements were 
made at Roosevelt Roads. The instrumentation and setup were similar to that described 
by Maring et al., [2000], however here we will only describe them briefly. 

A sample intake pipe brought air from the roof of the trailer laboratory down to an 
instrument table. Just above the table, ihe pipe made a 90 degree bend and ran 
horizontally to provide sample air to the instruments. At certain points, smaller tubes 
were connected such that the original air sample was teed off to different instruments. 

This intake quantitatively samples aerosols (less than 10% error) for up to 10 pm 
diameter (5 pm radius - hereafter, we will denote size in terms of radius), but has up to . . 
50% error for larger particles. The lowest meter of the tube includes an inline heater to 
control the relative humidity of the sample. The total flow rate in the intake was set to 90 
L. min' 1 , which provided a laminar flow such that diffusional losses for 3 nm particles 
were less than 4%. 

Aerosol number size distributions were measured with a series of instruments. A TSI® 
Aerodynamic Particle Sizer (APS) Model APS33 measured large aerosols from 0.4 to 
>7.5 pm aerodynamic radius with a resolution of 32 channels per half decade. Smaller 
particles (0.075 to 0.425 pm) were measured by a TSI® Scanning Mobility Particle Sizer 
(SMPS) Model 3934L. Corrections to the aerosol size distributions, for aerodynamic and 
diffusional loss, are discussed in Maring et al., [2000]. The combined information from 
these instruments give aerosol number size distribution in 150 radius bins ranging from 
about 0.007 pm to 7 pm, during a twenty minute measurement. From the size 
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distribution, we can calculate a particle effective radius (equation 5) in order to compare 
to MODIS. 

Aerosol light scattering measurements were performed using a TSI® integrating 
Nephelometer Model 3563 [Anderson et al., 1996]. The instrument measures total and 
backscattering by aerosols at 0.45, 0.55 and 0.70 pm. Data were corrected for angular 
non-idealities according to Anderson and Ogren [1996J. Data were recorded as five 
minute averages during the period 2 July to 24 July. 

Nearly continuous measurements by both sets of instruments were made during PRIDE. 
The data we use from the Particle Sizers were from the twenty minute measurements 
taken between ±1 hour of MODIS overpass, daily for the period 3 July to 14 July and 16 
July to 18 July (there were instrument problems on 15 July and after 18 July). The 
nephelometer data used also correspond to the same overpass time windows (but 
including 15 July and 18-24 July). 

4. Validation of MODIS Retrievals 

By deploying multiple instruments at different locations, we had a better chance of 
catching a comparison to MODIS. Also, more comprehensive studies of aerosol 
properties could be performed. Because data from each validating instrument were taken 
at different frequencies, slightly different approaches are used to compare MODIS with 
each data set. 

To compare MODIS optical depth with ground sunphotometer data, we use the spatio- 
temporal approach outlined in lchoku et al., [2002a], where the main assumption is that 
spatial statistics from satellite retrievals can be compared to temporal statistics from point 
observations. According to in-situ data taken at Roosevelt Roads, the average wind 
speed within ±1 hour of satellite overpass was similar to the 6.0 m/s used in the lookup 
tables (average = 5.98 m/s, standard deviation = 0.82 m/s). If we assume that the parcel 
can travel in any direction, aerosol can travel up to 22 km in any direction within an hour. 
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For easier computations, we chose a 50 km by 50 km grid as our MODIS validation box. 
That is, each validation box is composed of five 10 km x 10 km MODIS aerosol retrieval 
pixels in each direction (25 total). 

For the sunphotometer direct sun measurement data, we counted all valid observations 
within ±30 minutes of MODIS overpass. For both MODIS and ground-based 
measurements, we computed statistics for comparable quantities. Acceptable 
comparisons were defined when aerosol was retrieved in at least 20% of the pixels (i.e. 
five 10km x 10km pixels out of a possible 25) and at least two valid sunphotometer 
retrievals (out of a possible four or five for AERONET, and more for other types). 

When comparing size distribution retrievals, we were less stringent. Because only one 
almucantar was performed per hour, we could not usually expect to get two almucantars 
within the one hour period. Therefore, we relaxed to allow a single valid AERONET 
retrieval within a two hour period (±1 hour of overpass). For comparisons with in situ 
scattering and aerosol size parameters, we also used averages taken about ±i hour of 
MODIS overpass. However, we still used the 50 km x 50 km box for MODIS. 

4.1. Aerosol Optical Depth 

Total aerosol column optical depth at a particular wavelength, x x , is the quantity most 
easily “validated.” Figure 5 shows optical depth of MODIS compared to sunphotometer 
for two wavelengths (0.87 and 0.55 pm), separated by sunphotometer type (AERONET, 
MICROTOPS and AATS). Only at the 0.87 pm channel, can MODIS be directly 
compared with all three sunphotometers. At 0.55 pm, linear interpolation in log/log 
space was performed (between 0.50 and 0.67 pm for AERONET, between 0.440 and 
0.675 pm for Microtops and between 0.525 and 0.864 pm for the AATS). The lines are 
linear regression fits. We see that, for a given wavelength, there are no systematic 
differences between fits for different sunphotometers, implying that all sunphotometer 
data can be grouped together. 
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Now combining the data from all three sunphotometers into one data set, we co-locate 
and compare aerosol optical depth at 0.87, 0.66 and 0.47 pm (Figure 6). For each point, 
standard deviation for MODIS is plotted vertically, and represents the spatial standard of 
optica! depth within the 50 x 50 km box. Standard deviation for the sunphotometers, 
representing the temporal standard deviation within the one hour period, is plotted 
horizontally. We see that in general, the standard deviations are comparable in 
magnitude. Regression lines are given for each wavelength, and we notice that for all 
three wavelengths, the magnitude of the y-intercept is less than 0.04, implying little or no 
surface contamination to the retrievals. Both the 0.87 and 0.66 pm regression lines lie 
within the expected retrieval error over ocean (Ax = ±0.03 ± 0.05x - thin dotted lines), 
defined by Remer et al., [2002]. However, the 0.47 pm line lies outside (over-prediction 
by MODIS). Over all, we see that from long to short wavelengths (0.87pm to 0.66pm to 
0.47pm), MODIS goes from under-predicting to over-predicting the AOD. 

4.2. Angstrom Exponent and Spectral Dependence of AOT 

A parameter used to analyze the dependence of the optical depth on wavelength is the 
Angstrom exponent, a, defined as: 

ln(r A / t, ) 

(7) 

where x is the optical depth at wavelengths X, and [Eck et al., 1999]. Figure 7 plots the 
Angstrom exponent calculated from 0.66 and 0.47 pm for MODIS versus that for the 
sunphotometers (using the interpolated 0.47pm optical depth). Only the Angstrom 
exponent for optical thickness (at 0.66 pm) greater than 0.15 are plotted, because at low 
optica] thickness a small optical depth error at one or both wavelengths can introduce 
huge errors to the Angstrom exponent (as seen by Ignatov et al. [1998]). Yet, MODIS 
consistently over-predicts the Angstrom exponent. These results come from the spectral 
dependence of MODIS AOD retrieval quality. Also plotted in Figure 7, is a comparison 
of MODIS compared with the in-situ scattering estimates of Angstrom exponent. While 
this is only a single point (July 16), and it is at slightly different wavelengths, it is 


17 



consistent compared to the sunphotometer retrievals. All MODIS Angstrom exponents 
lie above the one to one line. 

Another way of looking at the spectral dependence of AOD, is to compare MODIS 
optical depth retrieval to sun photometer retrieval, plotted in wavelength space. Figure 8 
shows the spectral dependence of each retrieval for cases with x 670 greater than 0.15. For 
the same color curves, the solid curve connects the sunphotometer data, while the dotted 
curve connects MODIS AOD. On two dates, coincident sunphotometer measurements 
were made during MODIS overpass. On July 16, Roosevelt Roads (AERONET) and 
Roosey Roads (MICROTOPS) were operated next to each other, and observed nearly 
identical spectral optical thickness (red curves). On July 4, the AATS flew directly above 
the ship borne MICROTOPS, and observed similar optical depths in the red and near IR 
wavelengths, but varied somewhat at other wavelengths (dark green curves). For each of 
these pairs, there is only a single MODIS retrieval. The sunphotometer curves display 
very little slope in contrast to the MODIS retrievals. 

4.3, Aerosol Size Parameters 

The Angstrom exponent, a measure of the spectral dependence of aerosol optical depth, is 
related to the size distribution of the measured particles. Small Angstrom exponent 
values are associated with large particles in accordance with Mie scattering theory. The 
fact that MODIS over-predicts the wavelength dependence of optical depth, means that it 
probably under-predicts the size distribution. Figure 9 displays a comparison of the 
AERONET and MODIS effective radii (dots) calculated from Equation 5. Each point 
represents a successful sky radiance inversion from within two hours of MODIS 
overpass. When there were two retrievals (one before overpass, one after), they are 
marked with “A” and “B.” Note that there are fewer points than for direct optical depth 
comparisons, because of strict almucantar selection, and that size distribution retrievals 
are only possible from AERONET data. The actual optical depth (magnitude is 
represented by the “error bar” lines attached to each point) shows little influence upon the 
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quality of the retrieval. During PRIDE, all values of effective radius from MODIS are 
far lower than corresponding values from AERONET. 

As for the in-situ measurements of size at the surface, there are only two measurements 
that were coincident with MODIS overpass: 7 July and 16 July. Included on Figure 9 are 
the effective radii computed from the in situ size distribution (squares). Because surface 
dust measurements are related to total column measurements in the Caribbean [Smirnov 
et al., 2001a], we can plot the in situ measurements on the same graph with the total 
column measurements. Like AERONET effective radius, they are larger than MODIS, 
though somewhat closer. 

4.4. Errors of Validation 

Deeper analysis of the MODIS retrieval errors yields interesting relationships. We can 
summarize the previous plots by stating that the quality of optical depth retrieval is a 
function of wavelength, and that the spectral error gives rise (through Mie theory) to 
under-prediction of the particles’ size. Let us define an optical thickness error that is the 
departure from the expected errors (At = ±0.03 ± 0.05t) defined by Remer et al. [2002]. 
Thus, we are defining a measure that works for all values of optical depth, encompassing 
both the absolute error, 

absolute / \ 

Lrror — V^MODIS ‘ T SP ) 

and the relative error, 

relative , • \ l 

T error “ v T MODIS " T SP ) ' T SP > 

This “expectation error” can be defined as: 

expect _ (^MODIS ~ T SP ) 

error ( 0.03 + 0.05 r sp ) ’ 

where x SP is from a sunphotometer direct sun measurement. By defining this measure of 
the error, we can define a “large” expectation error when the magnitude is greater than 
unity, meaning that the actual error exceeded the expected uncertainty. For the effective 
radius, as we do not yet have an expectation of the error, so we will define the effective 
radius error simply as the relative error, that is: 
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RefCr = (Reff M0D ,s - Reff AERONET ) / Rcff AER0NE T > (9) 

where Reff AERONEX is retrieved from AERONET almucantar sky measurements [Dubovik 
and King, 2000]. For each wavelength, these quantities are plotted in Figure 10. This 
figure clearly illustrates what we learned from the spectral analysis. At 0.66 (xm, there is 
almost no relationship between t“^ c ' and Reff'^‘ :vc . The greater the over-estimation of 
the optical depth at 0.47 pm and under-estimation of the optical depth at 0.87 pm, the 
more severe the under-estimation of the effective radius is when compared to AERONET 
inversions. Even when all wavelengths of AOD retrieval fall within expectations, we still 
may have errors in effective radius. This is in contrast to the non-dust findings by Renter 
et al. [2002] that MODIS effective radii differed from AERONET values by only ±0.1 
pm. These wavelength dependence and effective radius retrieval problems seem to be 
unique to the dusty conditions during PRIDE. 

5. Discussion and Further study 

For PRIDE, we have compared MODIS aerosol retrievals over ocean with optical depth 
and aerosol size retrievals from sunphotometers. Regression lines of MODIS retrievals 
of optical depth (compared with sunphotometers) fall mostly within published estimates 
[Remer et al., 2002] at both 0.87 and 0.66 pm, but with 0.87 pm retrieved low and 0.66 
retrieved high. At 0.47 pm, MODIS significantly over-estimates the optical depth. These 
spectral discrepancies help to illustrate why MODIS retrieves small mode particles in 
dust regions, where we believe large particles dominate. These discrepancies were fully 
anticipated because of the necessity of assuming particle sphericity in the original 
MODIS algorithms. 

According to Dubovik et al. [2000], non-sphericity of dust particles can cause retrieval of 
a “spurious” small mode in the AERONET retrievals, especially when using radiance 
from scattering angles greater than 40 degrees. The actual dust phase function better 
resembles phase functions of smaller spherical particles than phase functions from 
spherical particles having dust size [Mischenko et al., 1997]. Figure 1 1 shows the size 
distributions retrieved from AERONET, for almucantars measured within two hours of 
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MODIS overpass. Nearly all of these AERONET retrievals show a mode centered 
around 0.1 pm, which is believed to be non-physical [Dubovik et al., 2000]. This brings 
up an interesting observation. That is, if AERONET has a spurious small mode, then the 
effective radius computed from AERONET is already too small, implying that the 
effective radii calculated from MODIS are actually worse. The in situ size measurements 
also suggest that MODIS retrievals may be far too small. The two squares plotted in 
Figure 9 are reasonably close to the one-to-one line. However, literature [Maring et al., 
2000] tells us that the Particle Sizers may under-sample the large particles (up to 50% for 
10 pm particles) so that in reality, the in situ effective radius would be larger. It is 
interesting to note that all three types of measurements seem to under-estimate particle 
size. 

According to Mischenko et al. f 1997], large scattering angles magnify optical differences 
between spheres and non-spheres. Therefore, for MODIS data, where scattering angles 
are all very large, we might expect to see some dependence of the retrieval upon 
scattering angle. Figure 12 shows errors of the spectral optical depth and reflectance 
fitting, against scattering angle (for x 660 > 0.15). Because the optical depth is large 
enough, we use the relative optical depth error, defined as 

('''MODIS ~ ' C AERONET V "''AERONET’ 

and the reflectance fitting error defined by Equation (4) in Section 2. The curves are fits 
to the data using second order polynomials. For all four wavelengths, the optical depth 
error is nearly constant up to a scattering angle of about 130 degrees. Then MODIS tends 
toward under-estimation of optical depth (especially at 0.87 pm) as the scattering angle 
increases. We also see that as the scattering angle increases, MODIS has more difficulty 
fitting the modeled reflectance to the observed reflectance. These trends, as functions of 
scattering angle, suggest problems with our aerosol model for dust-like aerosols. 

Problems could include our assumptions of size distribution and/or our assumptions 
about refractive indices. However, given the AERONET experience, errors are most 
likely due to assumptions of spherical particles, and the use of Mie theory in construction 
of the MODIS lookup tables. 
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PRIDE was our first opportunity to evaluate the performance of MODIS in dust settings. 
The PRIDE results tell us that the anticipated errors introduced by non-sphercity do not 
appear to significantly degrade our optical depth retrievals at 0.87 and 0.66 pm. 

However, non-sphericity does appear to affect optical thickness retrievals at other 
channels and causes severe under-prediction of the dust particle size. These problems 
can be corrected. Over time, the growing data base of co-located MODIS retrievals and 
sunphotometer measurements (in dust environments), will provide us with sufficient data 
to measure the phase function of the ambient, column integrated, non-spherical dust at 
MODIS observation angles. We plan to derive empirically corrected phase functions and 
to introduce these into the MODIS lookup tables. These new lookup tables should 
improve the ability of MODIS to observe and monitor dust aerosol. 
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Table la: Refractive Indices, Median, Standard Deviation and Effective Radius for 
small mode particles used in the MODIS Lookup Table 


a 

X=0.47-->0.86[im 
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Table lb: Refractive Indices, Median, Standard Deviation and Effective Radius for 
large mode particles used in the MODIS Lookup Table 
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Wet Sea salt type 

8 

1.53- 0.003i (0.47) 

1.53- 0.001i (0.55) 

1.53- O.OOOi (0.66) 

1.53- 0.000i (0.86) 

1.46-0.000i 

1. 46-0.001 i 

1.46-O.OOOi 

0.60 

0.60 

1.48 

Dust-like type 

9 



1.46-0.001i 

1.46-0.000i 

0.50 

0.80 

2.50 

Dust-like type 
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Table 2: MODIS overpass parameters for Roosevelt Roads, Puerto Rico, during 
PRIDE. 


Date 

DOY 

Time 

So!ar_Zen 

Sensor_Zen 

Relative_Azm 

Scat_Anqle 

Glint_Anqie 

Glint/NoGlint 

26-Jun-00 

178 

14:58 

20.89 

39.75 

25.11 

157.63 

59.21 

No Glint 

27-Jun-00 

179 

15:41 

12.08 

45.05 

148.31 

124.39 

35.23 

Glint 

28-Jun-00 

180 

14:46 

23.72 

54.76 

22.90 

146.19 

76.96 

No Glint 

29-Jun-00 

181 

15:29 

14.13 

25.18 

147.76 

142.16 

15.11 

Glint 

30-Jun-00 

182 

14:34 

26.59 

64.95 

20.98 

139.22 

90.00 

No Glint 

Ol-Jul-OO 

183 

15:16 

16.91 

2.05 

32.49 

164.78 

18.67 

Glint 

KEBHEDi 

184 

15:59 

7.92 

63.08 

129.94 

111.69 

58.18 

No Glint 


185 

15:04 

19.77 

29.38 

25.02 

166.02 

47.96 

No Glint 

04-Jul-00 

186 

15:47 

10.34 

52.09 

141.33 

119.60 

44.36 

No Glint 

05-Jul-00 

187 

14:52 

22.51 

47.56 

22.44 

152.20 

68.76 

No Glint 

06-Jul-00 

188 

15:34 

12.90 

36.37 

148.01 

132.27 

26.21 

Glint 

07-Jul-OO 

189 

14:39 

25.38 

60.22 

20.16 

142.94 

84.29 

No Glint 


190 

15:22 

15.64 

12.40 

153.34 

152.72 

7.14 

Glint 


191 

16:04 

7.79 

68.20 

123.00 

107.42 

64.12 

No Glint 

1 0-Jul-00 

192 

15:10 

18.50 

16.53 

24.11 

172.54 

34.24 

Glint 

n-jui-oo 

193 

15:52 

8.93 

58.50 

141.76 

114.34 

51.68 

No Glint 

lUPliiM 

194 

14:57 

21.37 

39.75 

20.92 

159.01 

60.11 

No Glint 

1 3-Jul-00 

195 

15:40 

11.56 

45.05 

150.25 

124.68 

35.38 

Glint 

1 4-Jul-00 

196 

14:45 

24.16 

54.75 

18.57 

147.50 

77.89 

No Glint 

1 5-Jul-00 

197 

15:28 

14.28 

25.18 

155.45 

141.40 

13.48 

Glint 

1 6-Jul-00 

198 

14:33 

27.01 

64.92 

16.63 

140.51 

90.94 

No Glint 

1 7-Jul-00 

199 

15:15 

17.10 

2.04 

25.28 

164.72 

18.96 

Glint 

1 8-Jul-00 

200 

15:58 

7.40 

63.05 

146.07 

110.74 

57.00 

No Glint 

1 9-Jul-OO 

201 

15:02 

20.02 

29.37 

18.50 

167.97 

48.73 

No Glint 

20-Jul-OO 

202 

1 5:46 

10.14 

52.08 

154.13 

118.69 

43.12 

No Glint 

2 1 -Jul-00 

203 

14:50 

22.94 

48.56 

15.97 

152.92 

70.82 

No Glint 

22-Jul-00 

204 

15:33 , 

12.91 

35.41 

159.70 

132.30 

23.67 

Glint 

23-Jul‘OO 

205 

14:38 

25.67 

60.19 

14.28 

144.32 

85.19 

No Glint 

24-Jul-00 

206 

15:21 

15.71 

11.49 

163.60 

153.08 

5.68 

Glint 

25-Jul-OO 

207 

1 6:03 

8.24 

67.50 

123.06 

107.85 

63.19 

No Glint 

26-Jul'OO 

208 

15:09 

18.66 

17.43 

14.96 

175.22 

35.77 

Glint 


DOY is the “Day Of Year” or the Julian date. The columns marked “Solar_Zen”, 
“SensorJZen”, Relative_Azm” and “Scat_Angle” denote the solar zenith angle, sensor 
zenith angle, relative azimuth (solar/sensor) azimuth angle and the scattering angle, 
respectively. 






































































































































































































































































Table 3: Identification of ground-based and airborne aerosol measurements 
coincident with MODIS overpass. 


1 MODIS 

AERONET Almucanturs 

MICROTOPS 

AATS 

In-situ 

Date 

Time (UTC) 

fflfl 

Location, Time (UTC) 

(Lat.Long), Location 

(Lat.Long) 


26-Jun-OO 

14:58 

129 





27-Jun-00 

15:41 

b 





28-Jun-00 

14:46 






29-Jun-00 

15:29 






30-Jun-OO 

14:34 

mm 

RR 14:27 




01 -Jul-00 

15:16 

mm 





02-Jul-00 

15:59 

mm 





03-Jul-00 

1 5:04 





Y 

04-Jul-00 

15:47 

n 


(17.961,-64.961), Chapman 

(17.93,-64.97) 

Y 

05-Jul-00 

14:52 

mm 


(18.167,-64.833), Chapman 


— 

06-Jul-00 

15:34 

mm 


(18.583,-65.167), Chapman 


BBS 

07-Jul-00 

14:39 

99 

RR 14:29, 15:29 

(18.250,-65.633), Roosy 


Y 

08-Jui-OO 

15:22 

B 




Y 

09-Jul-00 

16:04 

mm 

LP 14:35 

(18.250,-65.633), Roosy 


Y 

10- Jul-00 

15:10 

b 

LP 14:35 

(18.224,-65.646), Roosy 

(18.47,-65.15) 

Y 

11 -Jul-00 

15:52 



(18.224,-65.646), Roosy 

(18.56,-65.41) 

Y 

12-Jul-00 

14:57 

129 


(18.224,-6 5.646), Roosy 

(17.33,-65.43) 

Y 

1 3-Jul-00 

15:40 

mm 


(18.224,-65.646), Roosy 

(18.96,-65.13) 

Y 

1 4-Jul-00 

14:45 

MUM 




Y 

1 5-Jul-00 

15:28 

mm 


1— 


BH 

16 -Jul-00 

14:33 

mm 

RR 14:30, 15:30 

(18.224,-65.646), Roosy 


Y 

1 7-Jul-OO 

15:15 

mm 


(17.939,-67.017), One Day Cruise 


Y 

1 8-Jul-00 

15:58 

129 


(18.224,-65.646), Roosy 



1 9-Jul-00 

15:02 

9 

RR 14:30,15:30; LP 14:36 

(18.224,-65.646), Roosy 

(18.24,-65.47) 


20-Jul-00 

15:46 

■29 


i • . 

(17.58,-65.96) 


21 -Jul-00 

14:50 

■29 



(17.98,-65.65) 


22-Jul-OO 

15:33 

B 





23-Jul-OO 

14:38 

129 



(17.92,-65.38) 


24-Jul-00 

15:21 

zn 

■999H9H99I 



^Bl 

25-Jul-OO 

16:03 

mm 





26-Jul-OO 

1 5:09 

mm 






The column marked “G/NG” refers to whether Roosevelt Roads was within the MODIS 
glint mask. The listing for AERONET is for almucantars taken within two hours of over- 
pass. For the MICROTOPS, the location (and name of location) is exactly at overpass, 
while for the AATS-6, the location is the midpoint of its low-level leg. The last column 
denotes in-situ measurements within one hour of overpass. 
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Figure 1. MODIS “true-color” imagery seen on Jul 21, 2000. The image is a 
composite of Red Green and Blue MODIS channels. Note Puerto Rico to the left of 
center, and the glint reflection in the upper middle of the image. 


31 





Figure 2. MODIS combined land and ocean Optical Depth retrieval (at 550 nm) for 
July 21, 2000. Notice Puerto Rico to the left of the center. The glint region 
(reflection within 40 degrees of the specular reflection) is not retrieved. 
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Figure 3. Schematic showing how the field instruments were deployed together on 
July 4-6, where the filled circles, filled square and open triangle represent 
approximate locations of the AERONET, MICROTOPS and Ames sun- 
photometers. The solid line is a flight track of the Navajo aircraft, following the 
same heading as the MODIS satellite track (dotted). In situ data were taken at 
Roosevelt Roads. 
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Aerosol Optical Depth 


AERONET Optical Depth 
at 670 nm 

n.R 



Jul 1 Jul 15 Jul 30 


Day of Year (2000) 

Figure 4. Time series of columnar optical depth from Roosevelt Roads and La 
Paguera AERONET stations during PRIDE. Day 183 corresponds to July 1, 2000. 
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MODIS / Sunphotometer Optical Depth 

0.5 


0.4 


8 0.3 

< 

GO 

0 

O 0.2 
0.1 


0 

0 0.1 0.2 0.3 0.4 0.5 

Sunphotometer AOD 

Figure 5. Spatio-temporal comparison of MODIS optical depth retrievals (for 870 
nm and 550 nm) for each of the sun-photometers (AERONET, MICROTOPS and 
AATS) used during PRIDE. Lines represent linear regression by wavelength for 
each sun-photometer. The light gray line is the one-to-one line for comparison. 


Each Sunphotometer 



35 




MODIS / Sunphotometer Optical Depth 
All Sunphotometers 



Figure 6. Spatio-temporal comparison of MODIS optical depth retrievals for all 
sun-photometers grouped together. Solid, Dashed and Dotted are linear regressions 
for Blue, Red and IR wavelengths, respectively. The thin dotted lines are estimated 
errors (±0.03 ± 0.05r ) published by Remer et al. [2002]. 
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Angstrom Exponent 

Defined when Sunphotometer T g60 > 0.15 



Sunphotometer or InSitu Angstrom Exponent 


Figure 7. Comparison of MODIS Angstrom exponent to sunphotometer (dots) and 
to in-situ nephelometer measurements (square), showing that MODIS over-predicts 
for all cases (compared to the dotted one-to-one line). The text associated with each 
symbol represents the day of July (except 30 = June 30), location and type of 
instrument. Locations: R = Roosevelt Roads, L = La Paguera, O = Ocean; 
Instruments: C = Cimel (AERONET), M = MICROTOPS, A = AATS-6, 1 = In Situ. 
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Spectral dependence 
of aerosol optical depth 
where x „ > 0.15 

670 
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AATS Jul 21 
LaPaguera Jul 21 
Roosy Jul 1 6 
Roosevelt Jul 1 6 
AATS Jul 4 
Chapman Jul 4 
AATS Jul 23 

Roosevelt jun 30 

LaPaguera Jul 23 



0.1 W — Sunphotometer 
■— MODIS 


i i i 


400 600 800 1 000 

Wavelength (nm) 


Figure 8. Spectral dependence of PRIDE optical depth, when sunphotometer x m > 
0.15. For the same color curves, solid curves connect sun-photometer measurements, 
while dotted curves connect MODIS retrievals. Roosevelt Roads (AERONET) and 
Roosey Roads (MICROTOPS) wore operated next to each other on July 16, and the 
AATS flew above the RV Chapman (MICROTOPS) on July 4. Because each of 
these pairs of sunphotomer measurements occurred in the same MODIS “pixel,” 
only one MODIS retrieval is associated with each 
pai 
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Figure 9. Comparison of effective radius retrieved from MODIS versus AERONET 
(dots) and in-situ data at the surface (squares). The magnitude of each “error bar” 
is the optical depth (multiplied by 0.025). The numbers represent the date in July 
(e.g. 4 = July 4), and if followed by a letter, then it represents AERONET retrieval 
for the same date (e.g. “A” and “B” are two almucantars within one hour of MODIS 
overpass). The dotted line is one-to-one. 
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Errors of AOD and 
Effective Radius 



-100 -80 -60 -40 -20 0 

Effective Radius Error (%) 


Figure 10. Relationship between effective radius error and the optical depth 
expectation error (defined by Equation 9), for three wavelengths. Optical depth 
expectation error is defined so that a magnitude less than one (between thin dotted 
lines) represents a retrieval within the Remer et al. [2002] expected retrieval error 
(±0.03 ± 0.05 t ). A perfect retrieval would lie on the dotted zero line. The effective 
radius error is the relative error defined by Equation 9. 
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Volume 

Size Distribution 



Particle Radius (urn) 


Figure 11. Volume size distribution retrieved from AERONET almucantars (solid 
curves) and from in-situ data (dotted curves), corresponding to valid MODIS 
retrievals. For in-situ data, we are assuming an aerosol layer thickness of 1 km, to 
derive the same units. 
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Optical Depth Error (%) 


Scattering Angle Dependence 
of Optical Depth and Fitting Errors 
(where r > 0.1 5} 



110 120 130 140 150 160 170 


Scattering Angle (degrees) 

Figure 12. Optical Depth error (in percent) compared with scattering angle, for four 
MODIS wavelengths. MODIS reflectance fitting error is given by the Xs. The 
curves are second order polynomials 
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Popular Summary 

Mineral dust aerosols are airborne particles produced by wind erosion of the earth’s surface, and 
can be transported far distances by the prevailing winds. Their presence in the atmosphere 
perturbs global and regional climate, by scattering and absorbing sunlight. In large amounts, dust 
reduces visibility and adversely affects human health. A major source region for dust is the 
combined African deserts, from where summertime prevailing winds transport the aerosol 
towards the Caribbean Sea. In the summer of 2000 (June 26 to July 24, 2000), the Puerto Rico 
Dust Experiment (PRIDE) took place in Roosevelt Roads, Puerto Rico, to study the African dust 
aerosol transported into the region. 

Satellite systems have the ability to view large areas of the globe nearly simultaneously, thus are 
useful tools for monitoring large scale phenomena, including African dust transport. Recently, a 
new sensor, the MODerate Imaging Spectro-radiometer (MODIS), began retrieving global 
aerosol information, at much increased spatial resolution and spectral resolution (more 
wavelengths) from previous satellite missions. Over oceanic regions, this new capability allows 
for not only retrievals of the aerosol amount, but also determination of the aerosol type. PRIDE 
had the unique distinction of being the first major field experiment to allow for direct comparison 
of aerosol retrievals from MODIS with ground based and airborne based aerosol measurements. 

In this paper, we examine the ability of MODIS to retrieve dust aerosol properties over oceanic 
regions around Puerto Rico. The primary products from MODIS include the “spectral optical 
depths,” which are measures of the amount of light scattering and absorption in an atmospheric 
column, in turn related to the amount of aerosol, and the “effective radius,” a measure of the 
characteristic size of the particles, which is indicative of the aerosol type. Compared to ground 
based measurements, we show that MODIS accurately retrieves the aerosol optical depth at one 
wavelength, but systematically varies in quality when retrieving optical depths at other 
wavelengths. This discrepancy, related to wavelength, results in MODIS estimating an effective 
radius that is much smaller than the ground-based measurement. We believe that discrepancies in 
the spectral optical depth and dust size retrievals are, at least in part, associated with the use of 
only spherical particle models in the MODIS algorithm. Future generations of MODIS aerosol 
algorithms will include more accurate treatment of the non-spherical shapes of dust aerosol. 



